o 
o 



(N 



in 



c3 



fc-disjunctive cuts and a finite cutting 
plane algorithm for general mixed integer 

linear programs 



Markus Jorg 

Technische Universitat Miinchen, Zentrum Mathematik 
Boltzmannstra£e 3, 85747 Garching bei Miinchen, Germany 

joerg@ma.tum.de 



q ; July 26, 2007 

Abstract: In this paper we give a generalization of the well known split 
cuts of Cook, Kannan and Schrijver [5] to cuts which are based on multi-term 
disjunctions. They will be called /c-disjunctive cuts. The starting point is the 
question what kind of cuts is needed for a finite cutting plane algorithm for 
general mixed integer programs. We will deal with this question in detail and 
derive cutting planes based on ^-disjunctions related to a given cut vector. 
Finally we will show how a finite cutting plane algorithm can be established 
C ' using these cuts in combination with Gomory mixed integer cuts. 

O " 

o 

1 Introduction 

X 

In this paper we will deal with cutting planes and related algorithms for general mixed 
integer linear programs (MILP). As most of the results will be derived by geometric 
arguments we focus on programs that are given by inequality constraints, i.e. 



max cx + hy 

Ax + Gy <b 

xezp 



where the input data are the matrices A € Q mxp , G € Q mXlJ , the column vector b € 
and the row vectors c £ Q p , h £ Q g . Moreover we denote by the polyhedra P = {(x, y) : 
Ax+Gy < b} C W +q and P/ = conv ({(x, y) G P : x G W}) } C W +q the feasible domains 
of the LP relaxation and the (mixed) integer hull of a given MILP, respectively. We call 
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a MILP bounded, if the polyhedron P is bounded. We will also need the projection 
proj^ (P) := {x G MP : 3y G W 1 : (x, y) G P} of the polyhedron P on the space of the 
integer variables. 

By a cutting plane for P we understand an inequality ax + 0y < 7 with row vectors 
a G Q p , (3 G Q 9 which is valid for Pj but not for P. Using cutting planes gives a 
simple idea of how to solve a general MILP: Solve the LP relaxation of the MILP. If the 
optimal solution is feasible, i.e. satisfies the integrality constraint, an optimal solution is 
found. Otherwise find a valid cutting plane that cuts off the current solution and repeat. 
But unlike the pure integer case no finite exact cutting plane algorithm is known for 
general MILP. Therefor we remark that most cutting planes for general MILP such as 
e.g. Gomory mixed integer cuts [7] or mixed integer rounding cuts pj] are special cases 
of or equivalent to split cuts [5]. This fact and more detailed relations between these and 
other cuts are stated in [6]. Here a split cut is defined as a cutting plane ax + fly < 7 
for P with the additional property that there exists d G Z p , 5 G Z such that ax + (3y < 7 
is valid for all (x,y) G P which satisfy the split disjunction dx < 8 or dx > 8 + 1. So 
split cuts are defined not constructively but by a property, only. Now one can see in the 
following 'classical' example of Cook, Kannan and Schrijver [5] that split cuts are not 
sufficient for solving a general MILP in finite time. 

Example 1. The MILP 

maxy 

—x\ + y < 
-x 2 + y < 
xi + x 2 + y < 2 
x\, X2 G Z 

has the optimal objective function value but the problem cannot be solved by any algo- 
rithm that uses split cuts, only. A proof of this statement in a more general context is 
given in \Lemma 31 

On the other hand, as positive results in the context of cutting plane algorithms for MILP 
we can only give the following two special cases: For mixed 0-1 programs split cuts are 
sufficient for generating the integer hull Pj of a given polyhedron P. See e.g. in the 
context of mixed integer rounding cuts or [3] in the more recent representation of lift- 
and-project cuts. For general MILP, there only exists a finite approximation algorithm 
of Owen and Mehrotra [12] which finds a feasible e-optimal solution and uses simple split 
cuts, that means split cuts to disjunctions Xj < 5 V Xi > 8 + 1. 

So as split cuts fail in the design of a finite cutting plane algorithm for general MILP 
we want to generalize this approach to cuts that are based on multi-term disjunctions. 
Therefor we start in lsection 21 with the introduction of /c-disjunctive cuts and some of its 
basic properties. Afterwards we look at the approximation properties of the fc-disjunctive 
closures and deal with the question what kind of cuts is needed for an exact finite cutting 
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plane algorithm both in general and in special cases. Finally we derive a fc-disjunctive 
cut according to a given cut vector. In lsection 31 we turn to algorithmic aspects and give 
a way of how a finite cutting plane algorithm for general MILP can be designed using 
fe-disjunctive cuts in connection with the well known mixed integer Gomory cuts. Finally 
we will discuss the algorithm and give some interpretations. 

1.1 Preliminaries 

Here we repeat two basic results that we will need during this paper. The first one deals 
with the computation of the projection proj^ (P), the second one with the convergence 
of the mixed integer Gomory algorithm in a special case. 

Lemma 1. Let a polyhedron P = {(x,y) : Ax + Gy < b} be given. Then 

proj x (P) ={x £R P : v r Ax < v r b, Vr G R}, 
where {v r } re R is the set of extreme rays of the cone Q := {v € W 71 : G T v = 0,v > 0}. 

Proof. The statement follows by applying the Farkas Lemma, see e.g. |10j . 1.4.4. □ 

Next we look at the usual mixed integer Gomory algorithm [7]. Although the algorithm 
does in general not even converge to the optimum, the special case in which the optimal 
objective function value can be assumed to be integer, e.g. the case of h = 0, can be 
solved finitely using the algorithm. In detail we have the following 

Theorem 1. Let a bounded MILP JTJ) be given. Then the mixed integer Gomory al- 
gorithm terminates finitely with an optimal solution or detects infeasibility under the 
following conditions: 

1. One uses the lexicographic version of the simplex algorithm for solving the LP re- 
laxation. 

2. The optimal objective function value is integral. 

3. A least index rule is used for cut generation, i.e. the mixed integer Gomory cut 
according to the first variable Xj, that is fractional in the current LP solution, is 
added to the program. Here xq corresponds to the objective function value. 

Using the last theorem, it is obvious that we can check in finite time if there is a feasible 
point in a polytope with a given (rational) objective function value, as by scaling it can 
be always assumed that the optimal objective function value is integral. This is expressed 
in the following 

Corollary 1. Let a bounded MILP ([I]) with the additional constraint cx+hy = 7 be given. 
Then the mixed integer Gomory algorithm terminates finitely with a feasible solution or 
detects infeasibility. 
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(a) a 3-disjunction 



(b) a 4-disjunction 



Figure 1: Examples for /c-disjunctions in M? 

2 /c-disjunctive cuts 

2.1 Basic definitions and properties 

In analogy to the definition of a split cut based on a split disjunction we now define a 
/c-disjunctive cut that is based on a /c-disjunction that contains every integral vector. 

Definition 1. Let k > 2 be a natural number, d l ,...,d k £ TIP integral vectors and 
5 1 , . . . , S k £ Z. Then we call the inequalities d 1 x < S 1 , . . . , d k x < 5 k a k- disjunction, if 
for all x G V there is an i € {1, . . . , k} with d % x < 8 l . In this case we write D(k, d, 5) 
with d = (d 1 , . . . , d k ), 5 = (S 1 , . . . , 5 k ) for the k- disjunction. 

We note that we do not require the vectors d l ,5 l to be different. So every ^-disjunction 
is also a /c-disjunction for I < k. Especially every split disjunction is also a /c-disjunction. 
Moreover every fc-disjunction is a cover of Z p by definition. 

Definition 2. Let P C W +q be a polyhedron and ax + f3y < 7 be a cutting plane. Then 
ax + fly < j is called a k-disjunctive cut for P, if there exists a k-disjunction D(k, d, 5) 
with 

(x,y) G P : ax + (iy > 7 => d % x > <f , Vi € {1, . . . , k}. 

Of course every ^-disjunctive cut for P is valid for Pj by definition. According to the 
remark after IDefinition ll every ^-disjunctive cut is also a /c-disjunctive cut for I < k. So 
every split cut is a /c-disjunctive cut. 

Definition 3. Let P C W +q be a polyhedron. Then the intersection of all k-disjunctive 
inequalities is called the k-disjunctive closure of P and denoted by P^ . Analog the i-th 
k-disjunctive closure pj. of P is defined as the k-disjunctive closure of Pu ^. In the 
special case of k = 2 we will also write pW instead of plj, ■ 
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We want to remark that it is not evident if the /c-disjunctive closure P^ for a given 
polyhedron P is again a polyhedron in the case of k > 3. The both proofs of this 
property for the split closure [T] and [5] cannot be applied to the more general case. 
However, we will not further deal with this question, as our results in the following 
are independent of this property. We further remark that Definition 21 also applies in a 
natural way to closed convex sets P. Therewith it is guaranteed that the definition of 

(i) 

the i-th /c-disjunctive closure P^ of a polyhedron P is well defined. 

A valid cut to a given fc-disjunction can be computed as intersection cut to any basis 
solution of the LP-relaxation that is not contained in the disjunction according to [2]. 
In the case of k = 2, Andersen, Cornuejols and Li have shown [1] that intersection cuts 
are sufficient to describe all cuts to a given split disjunction. This result is not true for 
general /c-disjunctions. Here not every valid /c-disjunctive cut to a given disjunction is 
equal to or dominated by an intersection cut. This can be seen in the following 

Example 2. We look at the polyhedral cone C G M 2+1 with apex (3, g> o) ^at * s defined 
by 



-xx + y 


< 





-xi + y 


< 





xi+y 


< 


1 


X2+V 


< 


1 



Then y < is a 4- disjunctive cut for C to the 4- disjunction D := {x\ + X2>2,x\—X2> 
1, —X\ + a?2 > 1, —x\ — X2 > 0} . The set of all bases is given by any three of the above 
constraints. Computing the four relating intersection cuts to the 4- disjunction D we get 
that the point (g, 5 j g) * s valid for the four cuts and so y > g has to be satisfied. 

Although the properties of general ^-disjunctive cuts are more involved than in the case 
of split cuts, an investigation of these cuts is useful because every valid cutting plane for 
a given polyhedron is a /c-disjunctive cut for some k. 

Lemma 2. Let P C MP +g be a polyhedron and ax + (3y < 7 be a valid cutting plane. 
Then ax + (3y < 7 is a k-disjunctive cut for some k £ N. 

Proof. Let P be a polyhedron and ax + 0y < 7 be a valid cutting plane. The set M that 
is cut off by the above cutting plane is given by 

M := {(x, y) G P : ax + f3y > 7}. 

M contains no feasible points of Pi, So we have x 7L V for (x, y) G M. By lLemma li the 
projection of M can be expressed as 

proj x (M) = {x G W : A e x < b e , A l x < b 1 } 
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with wlog integral matrices A e , A 1 with rows dj satisfying gcd(aj) = 1 and integer vectors 
b e ,b l . We modify the coefficients of the vectors b e , b l by 

bj = [bf\ + l, 

b\ = m. 

Altogether we get that ax + (3y < 7 is a ^-disjunctive cut according to 
D(k,-(A e ,A l ),-(b e ,b 1 )). □ 



It is our goal to compute the mixed integer hull of a given polyhedron using fc-disjunctive 
cuts. Of course this should be done 'as simple as possible', what means that both the 
number k of hyperplanes needed for the disjunctions and the number of iterations in 
a cutting plane procedure should be small. At least the latter property can be easily 
realized as the next theorem shows. 

Theorem 2. Let P C W +q be a polyhedron. Then Pi = P$ . 



Proof. We will show that every valid inequality ax + (3y < 7 for Pj is a 2 p -disjunctive 
cut for P. This is sufficient for the theorem. Using [Lemma 21 we get that ax + f3y < 7 
is a ^-disjunctive cut with a related disjunction D(k,d,5). So the claim is shown for 
k <2 P . Otherwise the number of inequalities of the disjunction can be reduced until the 
required limit of 2 P . Since D(k,d,5) is a fc-disjunction we have 

Vx G Z p 3i G {l,...,fc} : d*x < S l . 

On the other hand we take the set of all integral vectors with the property d l x = 5 l 
for a given i £ {1, ••• , k}. Either it exists now a vector x G W with d l x = 5 l and 
d?x > 5i , Vj G {1, . . . , k} \ {i}, or we can expand the disjunction by setting the right 
hand side of the inequality to 5 l — 1 and repeat this consideration. This may also lead 
to the case that the inequality can be dropped. Therewith we can restrict ourselves to 
disjunctions with the additional condition: 

Vi G {1, . . . , k} 3x* G IP : d i x i = 5 1 A d j x l > 8 j , Vj G {1, . . . ,k}\ {i} 

The set conv ({x 1 , . . . , x k }) contains except for its vertices {x 1 , . . . , x k } no more integral 
vectors: Assumed that there was another integral vector z G conv ({x 1 , . . . ,x k }) we had 
d l z < 5 l for an i. This contradicts the definition of the vectors x l . So we have constructed 
a set that contains exactly k integer points as its vertices. This will lead to a contradiction 
for k > 2 P . Then we had at least two vertices v, w with the additional property that each 
component Vi,Wi, i G {1, . . . , n} of both vectors is either even or odd. So |(u + w) is an 
integral vector which is contained in conv ({x\...,x k }). This is a contradiction to the 
properties of the set. □ 



We look at an easy example to see that 2 p -disjunctive cuts are needed in general to 
compute the mixed integer hull of a polyhedron in one step. 
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Example 3. We take the p- dimensional unit cube C = [0; l] p and define the polyhedron 
Q by 

Q = \x : ax < max ax, a G { — 1, lj p |. 

x£C 

Next we embed Q in the W +1 and define the polyhedron 

P = conv U X Q ^ xeQ. 

Of course it is Pi = C and the only valid k-disjunction for the cutting plane x p+ \ < is 
defined by the facets of Q itself. 

As ITheorem~2l shows, the mixed integer hull of a general polyhedron can be 'easily' 
generated with 2 p -disjunctive cuts in theory. Of course for practical issues the use of 
disjunctions with an exponential number of defining hyperplane is very expensive. So we 
will deal in the following with the second question we mentioned above, i.e. what kind of 
fe-disjunctive cuts we need at least in computing the mixed integer hull using a repeated 
application of k-disjunctive cuts. 



2.2 Approximation property of split cuts 

Before we further analyze which cuts we need to solve a MILP exactly, we will deal 
with the approximation properties of ^-disjunctive cuts. We repeat that already using 
split cuts is sufficient to approximate the optimal objective function value of any MILP 
arbitrarily exact. Therefor look at the series (7®)»eN of objective function values that is 
given by 

7 (i) := max{cx + hy\ (x, y) € P W } (2) 

for an arbitrary objective function cx + hy that is bounded over the polyhedron P. In 
detail we get the following 

Theorem 3. Let P C W +g be a polytope, cx + hy an objective function, 7* = max{cx + 

hy\ (x,y) € Pi} and 7M as defined in J2|). Then for all e > there is an iq € N with 

| 7 (io) _ 7 *| < e 

Proof. A proof of this statement in a slightly different form using a repeated variable 
disjunction can be found in the paper [12] of Owen and Mehrotra. Moreover the algorithm 
in this paper also gives a constructive proof. □ 

As we can approximate any optimal objective function value arbitrarily exact using split 
cuts, the use of general /c-disjunctive cuts becomes necessary for determining exact solu- 
tions, only. Moreover we want to remark that in practical applications already optimizing 
over the first split closure often gives a good approximation of the optimal objective func- 
tion value. This was in detail investigated by Balas and Saxena [4] for instances from 
the MIPLIB 3.0 and several other classes of structured MILP. 
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2.3 Solving MILP exactly 



We now get back to the question what kind of cuts is needed to solve a general MILP 
exactly. As we will see, this depends on the structure of the projection of the solution 
space on the x-space of integral variables. For example the important special case of 
the solution space being a vertex can be solved just using split cuts. However, we will 
see that in general the required number of disjunctive hyperplanes is exponential in the 
dimension of the integer space. We start with the case that the solution set contains 
relative interior integer points. 

Theorem 4. Let P C W +q be a polyhedron, cx+hy an over P bounded objective function 
and 7* = max{cx + hy\ (x, y) E Pi}. If 

relint (proj x ({(x, y) E Pi : cx + hy = 7*})) n TP j= 

then there is a k € N with max{cx + hy\ (x,y) E P^} = 7*. 

Proof. If maxjcx + hy\ (x,y) € P} = 7* there is nothing to show, so let max{cx + 
hy\ (x,y) € P} > 7*. That means especially that hit (projj^ (M)) n TP = where 
M := Pi (1 {(x,y) : cx + hy = 7*} denotes the solution set. Moreover let x* E 
relint (proj^ (M)) n TP. To proof the claim we have to show that cx + hy < 7* is a 
split cut for one of the polyhedra p( h \ k € N. 

Let Ajx + Giy < bi denote these inequalities in the representation of the mixed integer 
hull that constrain the set M. With [Theorem 31 we get 

lim max{a/ ii x + gny\ (x, y) G P^} = bn. (3) 

fc— >oo 

Moreover x* lies in the boundary of the projection projjf (M^) of each set := 
P^ n {(x, y) : cx + hy > 7*}. As M^ k+1 ^ C there exists an inequality px < tt that 
is valid for all of the sets proj^ {M^ k ') with the additional property 

PX = 7T, Vx S W°ix C^0> (4) 

as x* € relint (proj^ (M)). If we combine ([3j) and ((4]) we get as direct consequence that 
cx + hy < 7* is a split cut to the disjunction D(p, n) for some 

After we have seen that split cuts are even sufficient for solving an important class of 
MILP exactly, we turn to the general situation. The idea for finite convergence using 
fe-disjunctive cuts in general consists of the basic principle that there has to exist a k- 

disjunction D so that the interior of the projection proj^ (M) of the solution set is not 

(i) 

contained in D. If no appropriate fc-disjunction exists for all closures P fc then we cannot 
achieve a finite algorithm using ^-disjunctive cuts. On the other hand, if this condition is 
satisfied for each face of the solution set, finite convergence can be shown in the general 
case. 
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Theorem 5. Let P C W +g be a polyhedron, cx + hy an over P bounded objective 
function, 7* = max{cx + hy\ (x, y) G Pj} and M := Pj Pi{(x, y) : cx + hy = 7*}. // there 
exists for both M and all its faces f G F with relint (proj^ (/)) (7 V = a k-disjunction 
Df(k,d,5) with the property 



then there exists a n G N with maxjcx + hy\ (x,y) G Pc } = 7*. 

Proof. We prove the claim by induction over the dimension I of the solution set M. We 
start with 1 = 0. In this case k = 2 can always be chosen and the result is a special case 
of ITheorem 41 We assume now that the claim is true for I — 1,1 £ N. 
So let M be the solution set of m&x{cx + hy\ (x, y) G Pi} and dim(M) = I. Moreover let 
for k G N exist a /c-disjunction D(k,d,5) for M according to the assumption. We proof 
that cx + hy < 7* is a k-disjunctive cut to the disjunction D for one of the sets pj^ . 

Therefor we show that it exists a n G N such that (pj, n ' C\{(x,y) : cx + hy > 7*})nD = 0. 
With the disjunction D we define the polyhedron Q := {(x, y) : dx > 5 A cx + hy = 7*} . 
As cx + hy = 7* is a supporting hyperplane of Pi, there exists (x, y) with cx + hy > 7* 
and x $l D such that each of the inequalities c^x + /i/y < 7/ defined by (x, y) and a 
facet / of Q is valid for Pi. All inequalities cjx + /i^y < 7/ at most support M in an 
under dimensional face. So it follows either by induction hypothesis or by ITheorem~3l 

(n) 

that the inequalities cjx + hjy < 7^ are valid for some P^ . Herewith, the condition 
(P^ n {(x, y) : cx + hy > 7*}) fl D = is satisfied and the theorem is proven. □ 

We remark that it is necessary to involve all the faces of the solution set in the last 
theorem, as the following example shows. 

Example 4. We define the polyhedron P C R 3+1 through the vertices 



For the objective function vector (0, 0, 0, 1) we get that M is contained in a split disjunc- 
tion, whereas the cut according to the face x% = is a 3- disjunctive cut, only. 

We now deal with the question what cuts we need to solve a general MILP. Therefor 
we use special sets that can arise as solution sets of MILP to give a lower bound of 
the required number of disjunction terms. The idea is based on a generalization of 
|Example 1| 



x G relint (proj x (/)) 



x #D f 



(0,0, 0,0), (2, 0,0,0), (0,2, 0,0) 
(0,0, 1,0), (2, 0,1,0), (0,2, 1,0) 
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Lemma 3. Let P C W +q be a polyhedron, cx + hy an over P bounded objective function, 
7* = max{cx + hy\ (x,y) G Pj} > max{cx + hy\ (x,y) G P} and M := Pj n {(x,y) : 
cx + hy = 7*}. If projjf (M) C MP has k facets with each containing a relative interior 
integer point, than max{cx + hy\ (x,y) G P}?li} > 7*j V« G N. 

Proof. Let proj^ (M) be given with relative interior points {x 1 , . . . ,x k } G IP that are 
contained in pairwise different facets. Then we have x %3 := \{x % + x 3 ) G int (proj^ (M)) 
for all i ^ j. By presumption there exists (x t3 ,y lJ ) G P with cx*- 5 + > 7*. 
Moreover at least one of the points [x %3 , y' 13 ) is not cut off by an arbitrary k — 1-disjunctive 
cut. So the cut is valid for the set Q 13 := conv ((x lJ , y 13 ), Pi). As each cut can be classified 
by this property we get that PWj Q %J — ^k-f ^ ' 1S c l ear that Cli^j contains a point 
(x,y) with cx + hy > 7*. As -f^-i satisfies all presumptions and the solution set M does 
not change, the proof follows by induction. □ 



Therewith we can show now that we need cutting planes to an in the dimension p 
exponential number of disjunctive terms to solve a MILP in general. 

Theorem 6. Let P be a polyhedron, cx + hy be an over P bounded objective function 
with max{cx + hy : (x, y) G P} = 7*. Then in general 

max{cx + hy : (x,y) G P^l 1+1 } > 7*, Vra G N. 

Proof. Using [Lemma 31 it is sufficient to give an integer polytope Q C W with p = n + 1 
and at least 2 n + 2 facets that contains no interior integer point but in each facet a 
relative interior integer point. Therefor we define Q as the set of all (x,x n+ i) G M n+1 
with the property: 

ax — 7r(a)x n+ i < 1, a G {±l} n 
< x n+1 < 2 

with ir(a) := \{i G {1, . . . , n} : aj = 1}| — 1. We show that Q has the desired properties. 
Its vertices are contained in the hyperplane x n +i = or x n+ \ = 2. For x n +\ = the 
related polytope is the n-dimensional cross polytope. For x n+ \ = 2 the related polytope 
is generated by the vertices 1 + (n — l)ui. The last property follows from the fact that 
ax < 1 is active for a vector ±Uj if, and only if ax < 1 + 2n(a) is active for 1 ± (n — 
by definition of ir(a). So Q is integer. 

Let (z, 1) G Z n+1 be given. We take the side constraint ax < 1 + 7r(a) with = 1 <^=^ 
Zi > 0, i G {1, . . . , n} and get 

az = N > ^ Zj > |{i G {1, . . . ,n} : o, = 1}| = 1 + 7r(o). 

l<i<n Zi>0 

So (z, 1) is no interior point of Q. Moreover we can see that (z, 1) G Q for z G {0, l} n 
and that (z, 1) is a relative interior point of the facet ax — 7r(a)x n+ i < 1 for aj = 1 <^=^> 
Zi = 1, i G {1, . . . , n}. As and 1 are relative interior point of x ra+ i = and x n+ \ = 2, 
Q has all properties. □ 
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(a) x 3 = (b) x 3 = 1 (c) x 3 = 2 



Figure 2: The set (J constructed in lTheorem 61 for p = 3, projected on the a^-space 

So we have proven that in general at least 2 P_1 + 2-disjunctive cuts are required to solve 
a MILP exactly in a finite number of steps. We remark that we have an upper bound 
of 2 P as shown in ITheorem~2l With this result we can see that an exact cutting plane 
algorithm gets in general very expensive as a large number of disjunctive hyperplanes 
has to be computed. Moreover we have not yet discussed how to determine a cut to the 
related fc-disjunction. As intersection cuts according to basis relaxations do not generate 
strong cuts in general, this is an important issue for practical applications. On the other 
hand we have seen that a wide class of problems can even be solved using split cuts. 
Moreover the convergence properties of /c-disjunctive cuts depend on the structure of 
the projection of the polyhedron and the objective function. This fact suggests to use 
information of the projection in cutting plane algorithms. 



2.4 Computing ^-disjunctive cuts 



At the end of this section we want to give an alternative to compute strong valid k- 
disjunctive cuts. Unlike the usual generating of valid cuts for MILP, we need as additional 
input the vector (c h) to which we want to cut. Moreover we restrict ourselves on shifted 
polyhedral cones P = {{x, y) : Ax + Gy < b} with apex (x*,y*),x* TP and assume 
that the function cx + hy attains its unique maximum at (x*,y*) with value 7*. In this 
situation we can describe (x*,y*) as polyhedron given by the (over-determined) system 

Using [Lemma 11 the projection of the above system on the a>space - that is equal to x* 
- is given by 

proj x (P r ) = L € R p : v r ( ^ \ x < v r ( _M , Vr G r\ (6) 
with R being the set of extreme rays of the cone 

Q=\ve R m+1 : ( ^ } v = 0, v > j . (7) 
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As the cone Q is rational, we can assume that the extremal rays v r are elements of the 
additive group 

Gp^ := ju;GQ m+1 :w( ^ } GZ m+1 j. (8) 

Therewith we can use the above polyhedral description of proj x {Py* ) to define a valid 
fc-disjunction for P, that does not contain the apex (x*,y*). We do this by rounding up 
the right hand sides of the defining constraints of proj^ (Pj*) in (|6j). 

Lemma 4. Let P,(x* ,y*),(c /i),7* and proj^ (-By*) as defined above. Moreover define 
for r G R 

d r '■= v l,...,m^ ~~ v m+l c 

$ r : = [v r i,..., m b ~ <+i7*J + 1 

with v r = (^i ... m , ^m+i)- Then D(\R\, — d, — 5) is a valid \R\- disjunction for P that does 
not contain (x*,y*). 

Proof. By definition it is maxjcx + hy : (x, y) G P/} < 7*. So there is an e > such 
that cx + hy < 7* — e is valid but not optimal for Pj . Moreover the inequalities defining 
proj^ (P 7 *_ e ) are given by v\ m A — v r m+1 c with right hand sides v\ m b — v r m+l {^f* — e). 
The set proj^ (P 7 *_ e ) contains no integer points, so for all iGff there is a r G R with 

d r x > v[^ m b - <+i7*- 

Therewith is follows that the polyhedron {x G MP : dx < 5} contains no integer point 
in its interior. This is equivalent to the set D(\R\, —d, —5) being a valid |-R|-disjunction. 
Moreover, as the right hand side of each defining hyperplane of the projection has been 
enlarged by the definition of D it is obvious that (x*,y*) is not contained in the dis- 
junction. This proofs the lemma. We remark that the above definitions and the proof is 
similar to lLemma 21 □ 



So we have found a /c-disjunction that can be used to cut off the current LP solution 
(x*,y*). As we have mentioned at the beginning of this section we want to cut to the 
vector (c h). We can do this now using the right hand side 5 of the disjunction D. As 
for v r m+1 > the value of 5 r depends on the objective function value, we can compute 
the objective function value that corresponds to the value of 5 r that we have got by the 
rounding operation. The inequalities of D whose right hand sides S r are independent of 
the value of 7 can be omitted in this considerations. Taking the maximum of the related 
objective function values for all constraints gives us a valid cut to the vector (c h). As 
the disjunction does not contain x* , we can ensure that the current solution is cut off. 

Theorem 7. Let P,(x* ,y*),(c h) , 7* , proj^ (P 7 * ) and D(\R\, —d, —5) as defined above. 
Let for r G R with v r m+l > 

r 5 r -v r i...., m b 
7 : = — — t 

V m+1 

and 7 = max{7 r : r G R, v r m+l > 0}. Then cx + hy < 7 is valid for Pi and cx* + hy* > 7. 
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Proof. The validity of the inequality cx + hy < 7 follows directly using ILemma 4l as it 
is a disjunctive cut according to D(\R\, —d, —5) by definition. Equally it follows that 
cx* + hy* > 7. □ 

As we have finished the derivation of the /c-disjunctive cut, we want to add some remarks. 
Using the projection as fc-disjunction, we solve the problem how to find a suitable k- 
disjunction for cutting in general. This relates both to the selection of the number k 
and the selection of the defining hyperplanes of the disjunction. Moreover we have seen 
in the last subsections, that using information of the projection can be useful. On the 
other hand, the projection that we use corresponds to the predisposed cutting vector. So 
the selection of a suitable ^-disjunction is partially shifted to the selection of the cutting 
vector. Here it is i.e. open how to choose cutting vectors to get deep cuts in general. 
However, for solving a given MILP we will see in the next section that this approach 
leads to a finite algorithm if we use the objective function vector. 



3 Algorithm 

We now turn to an algorithmic application of the previous results and want to present 
an exact algorithm that solves a bounded MILP in finite time. It is based on a series 
of mixed integer Gomory cuts that is mixed with certain fc-disjunctive cuts which are 
required as discussed in Isubsection 2.31 The fc-disjunctive cuts we use here are similar 
to the ones we introduced in Isubsection 2.41 using the objective function as the vector 
to which we cut. As the assumptions that we have made there for the fc-disjunctive 
cuts are in general not satisfied, we have to do some modifications. So we will define 
fe-disjunctive cuts over general polyhedra P for an arbitrary cut vector (c h). We discuss 
the details of the generalization. It is clear that the equalities ©, ([7|, ((8|) also describe 
the projection for P being a general polyhedron and 7* being an arbitrary value of the 
objective function. Even the derivation of a valid fc-disjunction and a valid fc-disjunctive 
cut, respectively, is true, if the value 7* of the objective function cx + hy is not optimal 
for Pi. However, for the application in the algorithm we will define a slightly weaker 
version of the |-R|-disjunctive cut that does not always cut off the current LP solution, 
but can be used more general. We do this in the next 

Theorem 8. Let P be a polyhedron and 7* such that cx + hy < 7* is valid for Pi but 
not for P . Define the \R\- disjunction D(\R\, —d, —5) using the equalities J5|), j6|), j7|), 
JEj) with 

d T '■= v \,..., m A — v m+l c 

■= K,..., m & - *4+i7*l 

and let 7 = max{7 r : r £ R, v r m+1 > 0} analog to I TheorernPA Then cx + hy < 7 is a 
valid cutting plane for Pi. 
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Proof. By assumption proj^ (P-y*) cannot contain an integral point in its interior. So 
rounding up the right hand sides gives a valid |i?|-disjunction. Therefor cx + hy < 7 is 
a valid cutting plane analog to the proof of lTheorem 71 □ 



We go on with the single steps of the algorithm. We start with the usual mixed integer 
Gomory algorithm as long as we get either a feasible solution of the MILP or a solution 
of the LP relaxation that has a lower objective function value. This happens in finite 
time by |Corollary l| if we restrict ourselves to polytopes. If the objective function value 
has decreased we can apply [Theorem 8l and compute a valid |i?|-disjunctive cut using the 
objective function as vector to which we cut. Now we can apply the Gomory algorithm 
to the modified program again, until either a feasible solution is found or the objective 
function value decreases, and use lTheorem 81 again. In this way we get an algorithm that 
finitely terminates with an optimal solution to the given MILP or detects infeasibility. 



The formal algorithm is stated in Algorithm 1 



Theorem 9. Let a bounded MILP ([TJ) be given. Then Algorithm 1 either finds an optimal 
solution or detects infeasibility in a finite number of steps. 



Proof. The proof follows immediately with the following two facts: Every while loop (16) 
to (27) has only finite many iterations by |Corollary 1 as P is bounded by presumption. 



Similarly the outer while loop (14) to (31) has only finite many iterations as the possible 
number of different values 7 is finite. □ 



Before we further discuss the algorithm we will give two examples. We start with repeat- 



ing Example 1 



Example 5. Let again the MILP 

maxy 

— x\ + y < 
-x 2 + V < 
xi + x 2 + y < 2 
x±,X2 G Z 

with the optimal solution (1)1,1) of the LP relaxation be given. The mixed integer 
Gomory cuts according to x\ and X2 are given by —x\ + 2y < and —X2 + 2y < with 

decrt 

1 -1 )y = 0, y > 



the new LP solution (|,|,|). As the value of the objective function has decreased, we 
compute 7 as in \Theorem 81 The extremal rays of the cone {(1 1 
0} are the three vectors 

(1,0,0,1), (0,1,0,1), (0,0,1,1), 
so the projection proj^ (P 7 ) is given by 



-X\ 


< 





-7 


-X2 


< 





-7 


X\ + X2 


< 


2 


-7 
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Algorithm 1 Exact cutting plane algorithm 
1: Input: bounded MILP p 

2: Output: "optimal solution (x*,y*)" or "problem infeasible" if no solution exists; 
3: 

4: (x*,y*) := argmax {ex + hy : (x, y) € P}; 

5: 7* := max{cx + hy : (x,y) G P}; 

6: 

7: if P = then 

8: "problem infeasible"; break 
9: end if 

10: if x* G IP then 

11: "optimal solution (x*,y*)"; break 

12: end if 

13: 

14: while x* IP do 

15: 7 — 7*; 

16: while 7* = 7 do 

17: Compute Gomory cut ct^x + a 2 y < j3 to P, (x*,y*) by least index rule; 

18: P := P n {(x,y) : a l x + a 2 y < /?}; 

19: (ac*, y*) := argmax {er + /iy : (x, y) E P}; 

20: 7* := maxjex + hy : (x, y) G P}; 

21: if P = then 

22: "problem infeasible"; break 

23: end if 

24: if x* G ZP then 

25: "optimal solution (x*,y*)"; break 

26: end if 

27: end while 

28: 

29: Compute 7 = max{7 r : r G R,v r m+1 > 0} according to lTheorem 81 for P, (c fa), 7*; 

30: 7* = 7; 
31: end while 
32: 
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Inserting the current value 7 = | of the objective function and rounding gives 
max r=1)2i 3 7 r = max{0, 0,0} = 0. After applying the related cut y < we get as new LP 
solution the feasible point (2, 0, 0) and the algorithm stops with an optimal solution. 

Second we show how the algorithm works for the example of Owen and Mehrotra [12]. For 
this ILP the usual mixed integer Gomory algorithm does not converge to the optimum. 

Example 6. Let the ILP 

max x\ + X2 

8x1 + 12x 2 < 27 
8x1 + 3x 2 < 18 
xi,x 2 > 
xi, X2 € Z 

with the initial LP solution 1) be given. After applying the first possible cut to x\ the 
value of the objective function decreases and we can go to the second step of the algorithm. 
As we have an ILP it is proj^ (P 7 ) = P 7 with P 7 given by 

-xi - x 2 < -7 
8x1 + 12x 2 < 27 
8x1 + 3x 2 < 18 
xi,x 2 > 

By rounding we get finally the valid cut xi + X2 < 2 that relates to the optimal objective 
function value. 

The result of this example is typical for applying the algorithm on ILP. In this case we 
have to presume that all input data is integral and the k-disjunctive cut to the objective 
function reduces to the Chvdtal Gomory cut of the objective function vector. 

Concluding we want to discuss the algorithm. We have seen in the last example that for 
an ILP the ^-disjunctive cut reduces to an integer Gomory cut to the objective function. 
So the whole algorithm can be seen as a variant of the pure integer Gomory algorithm in 
this case. The crucial fact for finite convergence of the integer algorithm is the possibility 
to add cuts both to the objective function and to each variable if they are not integral. 
Using /c-disjunctive cuts to the objective function we have now the possibility to add cuts 
to the objective function in the case of MILP as well. Therewith we obtain a convergent 
algorithm in analogy to the integer case. 

Of course the complex part of the algorithm consists in computing the fe-disjunctive cut 
as the number \R\ of extreme rays v r of the cone Q grows exponentially. So an efficient 
algorithm for computing the extreme rays of the related cone is required. Moreover 
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we have to ensure that the computed rays satisfy the integrality constraints, i.e. are 
contained in the group Qp 9 . Therefor we can presume in practical applications the 
coefficients of the matrix A and the vector c to be integer. Then the integrality constraints 
are satisfied, if all of the extreme rays v r are integer. However, we will not further deal 
with this issue here, but refer to the papers of Henk and Weismantel [9J and of Hemmecke 
[8] and the references therein. They state several algorithms for this and the similar 
problem of computing Hilbert bases of polyhedral cones. 

At last we want to give a further interpretation of the algorithm. Therefor we assume 
that the feasible domain P is full dimensional and bounded. One can see that in this case 
we can always choose an optimal solution of the MILP such that q defining inequalities 
of P are active. So the solution is contained in a (p + q) — q = p-dimensional face of P. 
Therefor we can solve the MILP by solving each of the related p-dimensional subproblems 
and taking the best solution. Moreover the set of feasible solutions in each p-dimensional 
face is discrete in general, so solving a MILP for a p-dimensional face can be interpreted 
as solving an ILP, as we could apply a suitable affine transformation. This means that 
solving a MILP can be seen as parallel solving of several ILP. Especially every valid 
cutting plane for Pj is even valid for each of the discrete subproblems. Therefor we need 
information of the related discrete subproblems if we want to generate strong valid cuts. 
As the number of p-dimensional faces of P grows exponentially, this interpretation also 
gives another reasoning that we need ^-disjunctive cuts with an exponential number of 
defining disjunctive hyperplanes to solve general MILP. Within the algorithm we can 
find the p-dimensional subproblems in the facets of the polyhedron proj x (Pj*), where 
the value of the right hand side S 1 can be related to the current objective function value 
of the subproblem. 
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